N to A electromagnetic transition form factors from Lattice QCD 
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The magnetic dipole, the electric quadrupole and the Coulomb quadrupole amplitudes for the 
transition 7A'^ A are evaluated both in quenched lattice QCD at /3 = 6.0 and using two dynamical 
Wilson fermions simulated at P — 5.6. The dipole transition form factor is accurately determined at 
several values of momentum transfer. On the lattices studied in this work, the electric quadrupole 
amplitude is found to be non-zero yielding a negative value for the ratio, Rem, of electric quadrupole 
to magnetic dipole amplitudes at three values of momentum transfer. 
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I. INTRODUCTION 

Recent photoproduction experiments on the nuclcon 
at Bates [J] and Jefferson Lab ^ have produced accu- 
rate measurements on the ratios of electric and Coulomb 
quadrupole amplitude to the magnetic dipole amplitude. 
Non-vanishing values for these ratios are thought to be 
connected with nucleon deformation. 

Deformation is a common phenomenon in nuclear and 
atomic physics. Classically, the multiphoton coincidence 
experiment of taking a flash photograph or observing 
an illuminated object distinguishes a deformed dumbcU 
from a spherically symmetric sphere. Quantum mechan- 
ically, a multiphoton coincidence experiment could also 
determine that a J=0 ground state of a diatomic molecule 
has a deformed shape. However, usually in electromag- 
netic probes of microscopic systems, we are constrained 
to make measurements associated with one-photon ex- 
change, corresponding to a matrix element of a one-body 
operator. In the case of a diatomic molecule, the one- 
body charge density of the J=0 state is spherically sym- 
metric, and cannot reveal the deformation that is present 
in the system. 

In many cases, however, when a nuclear or atomic sys- 
tem is well approximated by a deformed intrinsic state, 
it is still possible to observe its deformation using a one- 
body electromagnetic operator. We consider here the 
lowest order electric multipole, the quadrupole moment. 
For an axially deformed object, the quadrupole moment 
in the body-fixed intrinsic frame is given by 
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quadrupole moment is negative. For collective rotation of 
the deformed intrinsic state Q , the relation between the 
spectroscopic quadrupole moment, Q, measured in the 
laboratory frame and the intrinsic quadrupole moment, 
Qo, in the body-fixed intrinsic frame is given by 



- J{J+l) 
(J+l)(2J-h3) 



Qo 
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where p{r) is the charge density distribution. If Qq is 
positive, the object is prolate with the polar axis longer 
than the eauatorial axis. In contrast, for an oblate obiect 



where J is the total angular momentum of the system 
in the lab, K is the projection of J onto the z-axis of 
the body-fixed intrinsic frame, and we have considered 
the sub-state with azimuthal quantum number M=J. In 
the previous example of the J=0 diatomic molecule, al- 
though Qo 7^ 0, Eq. 2 yields Q=0 so that the deforma- 
tion, while present, is unobservable. Similarly, in the case 
of a nucleon with J=l/2, Q is zero although Qo may not 
be. However for the A with J=3/2, Eq. 2 shows that a 
deformed intrinsic state can be detected by the spectro- 
scopic quadrupole moment, Q. The E2 and C2 transition 
moments between the J=l/2 nucleon and J=3/2 nucleon 
have the same property of revealing the presence of defor- 
mation in the nucleon, the A, or both, and in this work 
we calculate these moments in lattice QCD and provide 
direct evidence for this deformation. 

The question whether the nucleon is deformed from a 
spherical shape was raised twenty years ago 4] and it 
is still unsettled. On the lattice, hadron wave functions 
obtained via density-density correlators can provide in- 
formation on the deformation of particles of spin higher 
than 1/2 H. This approach yields no information on the 
deformation of the nucleon for the same reason as the 
vanishing of its spectroscopic quadrupole moment. This 
is why in lattice studies, like in experiment, one looks 
for quadrupole strength in the jN A transition to 
extract information on the nucleon deformation. 

State-of-the-art lattice QCD calculations can vield 
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provide direct comparison with experiment. Spin-parity 
selection rules allow a magnetic dipole, Ml, an electric 
quadrupole, E2, and a Coulomb quadrupole, C2, ampli- 
tude. If both the nucleon and the A are spherical, then 
E2 and C2 are expected to be zero. Although Ml is 
indeed the dominant amplitude, there is mounting ex- 
perimental evidence over a range of momentum transfer 
that E2 and C2 are non-zero [ll Ql • A recent analysis of 
experimental results on the values of E2/M1 and C2/M1 
is shown to be incompatible with a spherical nucleon . 

Understanding the origin of a non-zero deformation 
is an important theoretical issue, which depends on 
QCD dynamics. The physical origin of non-zero E2 and 
C2 amplitudes is attributed to different mechanisms in 
the various models: In the constituent non-relativistic 
quark model the deformation was originally explained by 
the color-magnetic hyperfine inter-quark potential arising 
from one gluon exchange and producing a Z?-state admix- 
ture in the singlet-quark wave functions of the nucleon 
and the A The deformation due to the hyperfine in- 
teraction was also studied in 'relativized' quark models. 
In both cases the deformation obtained, measured by the 
ratio E2/M1, is smaller than that found experimentally. 

In the context of the constituent quark model, it was 
recently proposed that elimination of gluonic and quark- 
antiquark pairs leads to two-body contributions in the 
charge and vector current operators that produce a non- 
zero quadrupole moment in agreement with experiment 
using only s-wave functions for the nucleon and the A @ . 
In cloudy non-relativistic models @ as well as in chiral 
or cloudy bag models ^ the deformation arises because 
of the asymmetric pion cloud whereas in soliton models 
it is thought to be due to the non-linear pion field inter- 
actions. Most of the results obtained in cloudy baryon 
models predict values smaller than the experimentally 
measured ones suggesting that the deformation can not 
be entirely attributed to the pion cloud. This is agree- 
ment with lattice results where a non-zero deforma- 
tion is observed in the case of the rho meson even in the 
quenched approximation where pion cloud contributions 
are omitted. 

In the present work we will compare quenched and 
unquenched results for the transition matrix element 
^ A in order to examine sea quark contributions 
to the deformation 9]. An early, pioneering lattice QCD 
study 10] with a limited number of quenched configura- 
tions yielded an inconclusive result for the ratio of the 
electric quadrupole to magnetic dipole amplitudes, re- 
ferred to as EMR or Rem, since a zero value could not 
be statistically excluded. However the theoretical frame- 
work that it provided is still applicable and we will apply 
similar techniques to the present study making a number 
of improvements: (i) We use smearing techniques, which 
very effectively filter the ground state so that the time 
independent physical observables can be extracted from 
the correlators, (ii) The quenched calculation is done on 
two volumes at the same parameters to check the volume 



lattice allows us to simulate smaller quark masses, the 
lightest giving a ratio of pion to rho mass of 50 percent, 
(iv) For each lattice momentum transfer q we calculate 
the multipoles both in the rest frame of the nucleon and 
of the A. These two different choices of kinematics en- 
able the evaluation of the transition form factors at two 
different four-momentum transfers, (v) We study dy- 
namical quark effects by evaluating the form factors us- 
ing the SESAM configurations JJj that were produced at 
(3 = 5.6 on a lattice of size 16^ x 32 using Hybrid Monte 
Carlo for two degenerate flavours of dynamical Wilson 
fermions. vi) In all cases we use more configurations to 
improve the statistics. 

II. LATTICE MATRIX ELEMENTS 

The current matrix element for the 77V A transi- 
tion with on-shell nucleon and A states and real or vir- 
tual photons is shown schematically in Fig. ^ It has the 
form [13 

where p{s) and p\s') denote initial and final momenta 
(spins) and Ut{p',s') is a spin- vector in the Rarita- 
Schwinger formalism. 




A(p') N(p) 

FIG. 1: A''7 — > A matrix element. The photon couples to one 
of the quarks in the nucleon at a fixed time ti to produce a 
A. 

The operator can be decomposed in terms of the 
Sachs form factors as 

O^'^ ^ gMi{q^)Kl>;,+gE2[q^)Kl% + gc2{q^)Kl% , (4) 

where the magnetic dipole, Qmi, the electric quadrupole, 
Qe2, and the Coulomb quadrupole, Qc2, form factors de- 
pend on the momentum transfer = (p' — pY. The 
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and their expressions are given in ref. [IJ. The reason 
for using this parametrization for a lattice computation, 
as pointed out in ref. ^^Ii is that the Sachs form factors 
do not depend strongly on the difference between the 
nucleon and the A mass. From the Sachs form factors 
one can evaluate the ratio E2/M1 referred to above as 
EMR or Rem, and the analogous ratio of the Coulomb 
quadrupole amplitude to the magnetic dipole amplitude 
C2/M1, known as CMR or Rsm, which are the tar- 
get of recent experiments. Using the relations given in 
refs. 1301 the ratios Rem and Rsm in the rest frame 
of the A are obtained from the Sachs form factors via 



and 



Rsm — — 



|q| Gc2{q') 



2mA GMiiq'^) 
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The lattice construction of the appropriate matrix el- 
ements for the evaluation of these form factors follows 
closely that of ref. ^3|- The computationally most de- 
manding part in this evaluation is the calculation of the 
three-point correlation function that corresponds to the 
diagram shown in Fig. ^ and it is given by 



(G^"^(t2,ii;p',p;r))= expHp'-X2)exp(-fz(p'-p).Xi)r^"(f7|T[x?(x2,i2)j''(xi,ii)x''(0,0)] | ) , 

X2, Xi 

, (7) 

where an initial state with the quantum numbers of the nucleon is created at time zero and the final state with the 
quantum numbers of the A is annihilated at a later time ^2- The photon couples to one of the quarks in the nucleon at 
an intermediate time ti producing a A. The corresponding three-point correlation function for the transition A — > jN 
is given by 



{G^^'^ih,h;p',p;T))^ J2 cxpHp ' • X2) cxp(+*(p ' - p) • xi) r^"( 1] | T [x"(x2, t2)j^(xi, ti)xf (0, 0)] IQ] 

X2, Xi 

For the spin-i source , x''(x, i), and the spin-| source, x^^i^jt), we use the interpolating fields 

X^^x) = ^£"''{2 [u^'^ix] C^ad'ix)] u%x) + [u^'^ix) Cj^u\x)] d%x)] 
and for the projection matrices for the Dirac indices 



(8) 
(9) 
(10) 

(11) 



For large Euclidean time separations ^2 2> 1 and ti ^ 1, the time dependence and field normalization constants 
are canceled in the following ratio ^3 
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(G^^"^ft2, h;p ^ p; F)) {G^^'' ^{h,h; -p, -p Ft)) 
(5yG,0^(t2,p';F4)) (G^^(t2,-P;r4)) 



t2-tl>l,tl>l 



1/2 



n,(p',p;r;A*) 



where G and Gfj are the nucleon and A two point functions given respectively by 

(G~^(t,p;r)) = ^e-^P-F''" (f]|Tx"(x,t)x''(O,0) 

X 

(G^,^(t,p;F)) = ^e-*P-F^"(17|Tx^(x,t)x?(O,0)|17) 



(12) 
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The phase in Eq. E|is the same as that of G^-'^^(t2, ti', P p) since formany we have 



n,(p',p;r;Ai) 



Ea + TOA 



-1/2 



1 



^EN + ruNV'^' (G^^"^(t2,ti;p',p; r)) 



Ea J \ imj^J \ 2En 
We use the lattice conserved electromagnetic current, j^{x), given by 



(14) 



/ 



(15) 



symmetrized on site x by taking j'^ix) — > [j''(a;) + — /i)] /2, where Q/ is the charge of a quark of flavour / and 
K/ is its hopping parameter. 

In the nucleon laboratory frame p = and p ' = q. The Sachs form factors can be extracted from the plateau 
values of no-(p ', p ; F; fj,) for specific combinations of matrices T and A vector indices a. The expressions for general 
momentum transfer q are obtained using the standard Euclidean non-relativistic representation for the 7 matrices [T^ 
with ei234 _ rpj^j, kincmatical functions K'^^ in Euclidean space are given by 



"■Ml ^ 



{mA + itin)^ + 2mN 



(mA+mjv) ^^r^.c.l3py, 



2mAr 

2m AT 



2mAr 



(16) 



with rJ(Q^) = [(toa + mjy)^ + Q^] [{niA — rnM^ + ^"^^ Q = q^ Q"* = is the lattice momentum transfer giving 
= —q^. By and p' ^ we now denote Euclidean space momenta defined analogously to Q^. The Rarita-Schwinger 
spin sum for the A in Euclidean space is given by 



^Ua{p,s)Ur{p, S) 

s 

and the Dirac spin sum 



-i"f.p + TOA 
2mA 



^ u{p, s)u{p, 



2p„Pr .Pair - Prla 1 
5a T + -17-^ I 5 ^7<t7t 

dm^ dTOA o 



~i"f.p + TOjv 



2mN 



(17) 



(18) 



We generalize in what follows the expressions of ref. 10] to allow momentum transfers in any spatial direction. By 
selecting the time component of the current the Qc2 form factor is extracted from 



GC2 = C(q2) 



2 Toa 



SikqZ - qkqi{\^2EA/mA) 



C(q2) = 



3 AEatiin ( Ea + m-A 



2 toat + TOA V Ea 



1/2 



nz(q,0;zrfe;4) 

^2 \ 1/2 



3to^ 



(19) 



where Ea = \/p' ^ + to^ and the indices fc and I denote spatial directions. By selecting the spatial component of the 
current the Qm\ and Qei form factors are extracted from 



Qui = C(q2)e-"=4j_n^(q^o.^p^.;) 



q'^Wi (q, ; r, ; /) - g'Hfc (q, ; r, ; fc) - ^ ( g'^Hfe (q, ; ; - q'^i (q, ; ; fc) 



and 



C(q2 



1 r_fcTT /_ 



Jtt /_ 



(20) 
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provided q'^ ^ (I ■ If we consider a momentum transfer that has zero component along the current direction Eas. I2UI 
and 12 II simplify to 



( Qui \ 


c(q2) r 


\ iGE2 ) 


qk 



n,(q,0;rfc;0 + ^nfc(q,0;r,;0 
Another possibility is to extract the transition form factors by using, instead of Ra, the ratio 

(G^^-^^(t2,ii;p',p; r)) 



i?W(t2,ii;p',p;r;/i) = 

or equivalently 

i?(,2)(i2,ii;p',p;r;Ai) = 



(<5„G^^(<2,P';r4)) 



(G^^(2ii,p;r4)) / (<5,,G,^^(2ti,p';r4)) 



1/2 



(G^^"^(t2,ii;p',p; r)) 



(G^^^(t2,P';r4)) 



(<5.,G^^(2ti,p;r4)) / (G^^^(2ii,p';r4)> 



1/2 



(22) 



(23) 



(24) 



In this work we choose the current along the z-direction 
and consider momentum transfers along the x-axis. In 
particular we consider the lowest allowed lattice momen- 
tum transfer q = (27r/A^a, 0, 0), where a is the lattice 
spacing and iV the spatial lattice size, with the exception 
of the large quenched lattice where we also consider mo- 
mentum transfer q = {Air/Na, 0, 0). For our choice of the 
momentum q it is reasonable to take in Eas. 1121 1231 and 
1^ the more symmetric combination 3/2{G^^ + G^'^) 
instead of SijG^^. With this replacement the second 
square root involving q^ in the overall factor G(q^) given 
in Eq. E|is absent. Using the fact that the momentum 
transfer has a component only in the x-direction Eas. 1191 
and|23simphfy to 



Ma) 



Co^^n4(q,0;-^r,;4) 

q"^ Ea 



(25) 



.(b) 
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Go 
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2rnA 
2rnA 



n2(q,0;zr2;4) 
n3(q,0;zr3;4) 



and 



^Mi = ^0 An2(q,0;ir4;3) 



(26) 



(27) 



Mb) _r 



n3(q,0;ri;3)-^ni(q,0;r3;3) 



(28) 



,(6) 
^E2 



Co 
3 



n3(q,0;ri;3) 



rriA 



ni(q,0;r3;3) 



(30) 

where Go is obtained from G(q^) given in Eq.^lby omit- 



lii 



ting the second square root [la| • When the A is produced 
at rest, the factors mA/i?A in Eas. 1251 1^ and QUI are ab- 



sent and Go ^ i/§ / i^iv+miv The resulting 

formulas then agree with those given in ref . . 
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FIG. 2: The ratio _R2(t2, ti; q, 0; ir4; 3). Data with a photon 
carrying momentum q and — q averaged over 50 quenched 
configurations are denoted by the stars. Data with a pho- 
ton carrying either positive or negative q averaged over 100 
configurations are shown by the x-symbols and the crosses re- 
spectively. These results are for a lattice of size 32^ X 64 at 
K — 0.1558. The current couples to the quark at ti/a = 8. 
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Smearing is essential for achieving ground state domi- 
nance before the signal from the time correlators is lost 
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FIG. 3: C/mi in units of natural magnetons using the ratios 
Ra (x's), R'^^^ (diamonds) and R^y'' (crosses) with the nucleon 
at rest (top) and with the A at rest (bottom) for dynamical 
quarks at k = 0.1565. The current couples to the quark at 
ti/a = 6. 
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FIG. 4: i'^ units of natural magnetons extracted from 

Eq. 1211 using the ratios R^ (x's) , r'^^^ (diamonds) and R''/^ 
(crosses) with the nucleon at rest (top) and with A at rest 
(bottom) for dynamical quarks at «; = 0.1565. The current 
couples to the quark at ti/a = 6. 
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Wuppertal smearing, d{x,t) — > c?™°^' (x, t), at the source 
and the sink. We smear the fermion interpolating fields 
according to 0, 

(x, i) = E ^' U{t))d{z, t) (31) 

z 

with the gauge invariant smearing function constructed 
from the hopping matrix H: 

F(x, z; U{t)) - (1 + ai/)"(x, z; U{t)), (32) 

where H{x,z;U{t)) - ELl(u^{^,t)S^.y-^ + U}{^ - 

j, i)(5x,y+i^ • The parameters a — 4.0 and n — are 

tuned such as to optimize the overlap with the baryon 
states. Quark propagators with a photon insertion are 
computed with the sequential source technique. There- 
fore for the three-point functions we require two inver- 
sions with the second inversion having a momentum de- 
pendent source. The sequential source technique requires 
that the photon couples to the quark at fixed time U 



ground states are identified. For the lattices used here we 
know from the nucleon and A two point functions that for 
ti > 5a the excited state contributions become negligi- 
ble. We also make use of the equal weighting of {C/j and 
{[/*} gauge configurations in the lattice action [l^l and 
the parity symmetry of our correlators to improve the 
plateau behaviour of the reduced ratio ncr(p',p ;r;/i). 
For implementation of this procedure we need an addi- 
tional sequential propagator with the photon carrying 
momentum — q. To see the improvement we compare in 
Fig. 121 the results obtained when for each configuration 
we consider the photon carrying either momentum q or 
— q to those obtained for the photon with momenta q 
and — q for equal statistics. As it can be seen, for large 
time separations, the quality of the plateau obtained for 
the reduced ratio Ila{p',p ;r;/x) is improved when the 
equal reweighting of {U} and {U*} is implemented en- 
abling us to fit over a larger range, far from the time 
insertion of the current. In all the results presented here 
this benefit is utilized: thus for each configuration the se- 
quential propagator is inverted twice, once with the pho- 
ton carrying momentum q and once carrying momentum 
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earlier plateau, which in this case is indispensable since 
the signal becomes too noisy beyond time separations 
ti — t2 > 6/a. Eg. 1301 provides an alternative way to ex- 
tract the electric quadrupole form factor. However the 
plateaus for the reduced ratios involved in the extraction 
of CJ^2 deteriorate when the nucleon is at rest leading to 
unreliable results. When the A is at rest, Eq. 1201 pro- 
duces a good plateau for the ratio Rc yielding results for 
consistent with those for Q^^^ ■ Therefore in what fol- 
lows, Eq. 1201 will only be used when the A is produced at 
rest to check consistency with the values obtained from 

Eq. [501 The ratio Rem is evaluated using ^^'^ ^i/i 
since this is applicable for both types of kinematics. Sys- 
tematic errors and noise become worse for Gc2- When 
the A is produced in motion the ratio Ra- cannot be used 
since the correlators G^^^ and G^^'^ do not agree even 
in the sign for time separations t2 — ti < 11 /a and af- 
ter that they are too noisy to be usable. This is seen 
in Fig. 101 where we show results for Gc2- ^^'^ other two 
combinations that yield Q^^^ show at best equally poor 
results. When the A is produced at rest the signs are 
consistent and the ratio Ra- , although noisy, can be used 
to look for a plateau. In what follows we will attempt to 
fit over the plateau range only for the case when the A is 
at rest in order to provide an estimate on Gc2 wherever 
possible. 



FIG. 5: in units of natural magnetons extracted from 

Eg. 1251 using the ratios (x's), R'^^ (diamonds) and R''^^ 
(crosses) with the nucleon at rest (top) and with A at rest 
(bottom) for dynamical quarks at k = 0.1565. The current 
couples to the quark at ti/a = 6. 



at large time separations from the source outweighs the 
additional costs. 

The ratios given in Eqs. El 1221 and provide three 
ways for extracting the form factors. In Fig. (21 we com- 
pare these three possibilities for evaluating Gm\- The 
plot illustrates the results for the unquenched case at 
K = 0.1565 - a similar behaviour is observed in all cases. 
We consider two kinematically different cases: one where 
the nucleon has zero momentum and therefore the A car- 
ries momentum q and the other where the A is produced 
at rest and therefore the nucleon has momentum — q. 
What is clearly seen for both kinematics, is that R^ as 
given in Eq.^l yields the best plateau, which starts as 
early as two time slices away from the time where the 
current couples to the quark. In contrast, the other two 
definitions require five time slices to show convergence 
to the same value. Evidently excited states contribu- 
tions come with the opposite sign in the matrix elements 
'yN A and A '^N canceling to a large extent in 
the ratio R^- The same conclusion is reached for 
extracted from Eq. [2H1 (without the toa /Ea factor when 
A is at rest). In Fig.^lwe show the analogous results for 
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FIG. 6; Gmi in units of natural magnetons (e/2mjv) for the 
unquenched theory on a lattice of size 16^ x 32 at /t = 0.1565 
and momentum transfer q = (27r/16a, 0, 0). Filled triangles 
and crosses denote results for g^^l {Eq^ and (Eal^ 
respectively when the nucleon at rest. The x's and stars de- 
note respectively results for G^j^j^ and G'^^ji when the A is pro- 
duced at rest. They are shifted to the left and right of the 
results obtained when the nucleon is at rest for clarity. The 
photon is injected at ti/a = 6 as shown by the arrow. 

In the case of the dominant amplitude Ml definitions 



,(a) 



and C/ltl can be used for both our kinematics and 



the results for Gm^i and G^'i are directly compared in 
Fig. 101 for the SESAM lattice at k = 0.1565 at momen- 
tum transfer q = (27r/16a, 0, 0) for both when the A is 
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FIG. 7: Ge2 in natural magnetons for the unquenched theory 
on a lattice of size 16'' x 32 at k = 0.1565 and momentum 
transfer q = (27r/16a, 0, 0). Filled triangles denote the results 
for when the nucleon is at rest. The x's and stars denote 
respectively results for t/^j 8-nd G^^l when the A is produced 
at rest. The rest of the notation is the same as that of Fig.|S| 
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t2 — ti > 4/a suggesting that contamination due to ex- 
cited states is more severe when the A is produced in 
motion. 

Figs. [3 and |S1 show the corresponding resuhs for Qe2 
and Qc2- For both form factors the data become very 
noisy for time separations t2 — ti > 7 /a. However for the 
case of the electric quadrupole we have a plateau region 
extending over four time slices enabling us to extract a 
value for Q^^2 both for A static and for A carrying a 
momentum. When the A is produced at rest the re- 
duced ratio needed for the extraction of Cy^2 shows simi- 
lar plateau behaviour as that obtained for tj^j whereas, 
as we already mentioned, when the A is not at rest this 
correlator is too noisy to be useful especially for the small 
quenched and unquenched lattices. The identification of 
the plateau region becomes particularly difficult for Qc2 ■ 
In Fig. |51 we show results for the ratio Ra only for the 
case when the A is produced at rest since, as we have 
already discussed, when the A has non-zero momentum 
the three point functions fiuctuate in sign making the 
ratio Ra unusable. From Fig. |S1 it can be seen that the 
results obtained for j involving matrix elements with 
the A vector index a in the same direction as that of 
the momentum transfer, have overall smaller statistical 
errors as compared to those for QqI- In fact, for most 
cases, the ratio Ra- from Eas. 1261 can only be determined 
for small time separations, t2 — ti, insufficient to neglect 
contributions from excited states. Therefore in what fol- 
lows we will mostly use Eq. 123*0 estimate Gc2- 

The behaviour of the transition form factors shown in 
Figs. El □ and El for the SESAM lattice at /t 0.1565 
is typical and it is observed for the other lattices and k 
values. 



FIG. 8: Gc2 in natural magnetons for the unquenched theory 
on a lattice of size 16"' x 32 at k = 0.1565 and momentum 
transfer q — (27r/16a, 0, 0). We show only results for the A at 
rest. The x's denote the results for (without the factor 
m^/EA in Eq. I25L stars and diamonds denote results for 
GqI using Eq. 1261 with the A vector index cr = 2 and a — 3 
respectively. The rest of the notation is the same as that of 
Fig.El 

the A rest frame Gm\ ^^'^ perfect agreement 

with a plateau region that sets in as early as two time 
slices from the insertion of the current. When the A car- 
ries momentum, G^^\ and agree for time separations 



III. RESULTS 

The quenched calculation of the transition form factors 
is carried out at (3 ~ 6.0 using lattices of size 16'^ x 32 
and 32'^ x 64. For the unquenched calculation we use 
the SESAM configurations a,t (3 — 5.6 on a lattice of 
size 16^ x 32. We use Wilson fermions with hopping 
parameters k given in TableQ]where we also list the values 
of the ratio of the pion mass to the rho mass. For all 
configurations, the number of which is given in Tabled 
we double the statistics performing the calculation both 
for q and for — q. 



To set the lattice spacing a in the quenched theory one 
can use the well known value of the string tension. How- 
ever since we want to compare quenched and unquenched 



both cases. Since we are calculating matrix elements in 
the baryon sector it is more appropriate to use the value 
extracted from the nucleon mass in the chiral limit to set 



TABLE I: K values and momentum transfers used for the 
evaluation of the transition form factors. The ratio of the rho 
mass to the pion mass at these k values is also given. We used 
the nucleon mass to set the scale. 





Quenched (3 — 6.0 16 


3 X 32 q 


2=0.64 GeV^ 


(GeV^) 




7TI-JY j TXlp 


m^(GcV) 


# of confs 


p at rest 


A at rest 










0.57 


0.64 


0.1530 




0.877(3) 


100 


0.55 


0.64 


0.1540 


0.78 


0.736(2) 


100 


0.50 


0.64 


0.1550 


0.70 


0.604(2) 


100 


0.40 


0.64 


«:c = 0.1571 








extrapolated 




Quenched /3 = 6.0 32^ 


X 64 


= 0.64 GeV^ 


(GeV") 


K 


Imp 


m^(GeV) 


# of confs 


p at rest 


A at rest 










0.50 


0.64 


0.1550 


0.69 


0.598(4) 


100 




Quenched /3 = 6.0 32^ 


X 64 q^ 


= 0.16 GeV^ 


(GeV") 


K 


mT,/mp 


m^(GeV) 


# of confs 


p at rest 


A at rest 










0.13 


0.16 


0.1554 


0.64 


0.537(4) 


100 


0.10 


0.15 


0.1558 


0.59 


0.469(4) 


100 


0.057 


0.13 


0.1562 


0.50 


0.392(4) 


100 


0.064 


0.13 


Kc = 0.1571 








extrapolated 




Unquenched (3 = 5.6 16^ x 32 q 


2 = 0.54 GeV^ 


(GeV") 


K 


ni^/mp 


m^(GeV) 


# of confs 


p at rest 


A at rest 










0.48 


0.54 


0.1560 


0.83 


0.837(6) 


196 


0.48 


0.54 


0.1565 


0.81 


0.755(10) 


200 


0.45 


0.54 


0.1570 


0.76 


0.655(9) 


201 


0.45 


0.54 


0.1575 


0.68 


0.527(8) 


200 


0.40 


0.53 


Ate = 0.1585 








extrapolated 



the nucleon mass is a^^ = 2.04(2) GeV (a = 0.089 fm) 
and in the unquenched a^^ — 1.88(7) GeV (a = 0.106 fm) 
with a systematic error of about 15% coming from 
the choice of chiral extrapolation ansatz. We note that 
if one uses the rho mass at the chiral hmit one obtains 
a-i=2.3 GeV(a = 0.087 fm) for both the quenched and 
the unquenched theory (IL ilft . .2Qi] , with a systematic er- 
ror of about 10% due to the choice of fitting range and 
chiral extrapolation ansatz . Therefore using the nu- 
cleon or the rho mass to set the scale yields results for 
a within the systematic uncertainties. We would like to 
stress here that in the ratios Rem and Rsm, which are 
the experimentally measured quantities, the lattice spac- 
ing does not explicitly enter and therefore a precise deter- 
mination of a is less crucial here than in other studies. In 
Table U we use the value for a obtained from the nucleon 
mass to convert the momentum transfer q at which the 
transition moments are evaluated to physical units. The 
lattice four- momentum transfer, Q^, depends on the k 
values through the nucleon and A masses. These values 
are given in Tabled As already mentioned we consider 
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FIG. 9: 



,(a) /^(a) 

(x-symbols) and 32^ x 64 (stars) at k 



in units of natural magnetons and 



CJ^^V^Mi in the quenched theory for lattices of size 16^ X 32 



0.155 and momentum 
transfer q = (27r/16a, 0, 0). For this comparison the data for 
the large lattice are shifted by two time slices to the left so 
that the coupling of the photon to the quark occurs at the 
same time point on the graph. The dashed and solid lines 
are the plateau values for the small and larger lattices respec- 
tively. 



at rest and the other with the nucleon at rest. For all k 
values, including Kc, for the case when the A is produced 
at rest the variation in is negligible. For the case 
when the nucleon is at rest the variation in over the 
range of n values considered here is sizable, in particular 
when q = {271 /'iia, 0, 0) where changes by a factor of 
two. We will come back to this point in section IV when 
we discuss the chiral extrapolation of the results. In all 
the evaluations we take u and d quarks of the same mass. 
The current couples to the quark at time slice ti = 6a, 
for the lattices of size IQ^ x 32 and at t\ = 8a for the 
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TABLE II: The plateau values of Gm\, GeI units of e/2mN and Rem in % for different fit ranges axe given together with 
the x^/d o.f. for the quenched and SESAM lattices in the rest frame of the A. Q^^l is also given for the cases where a plateau 
can be identified. With star we indicate the value that we adopted. 



16^ X 32 quenched 


16^ X 32 unquenched 


32^ X 64 quenched 


K = 0.153 


K 


= 0.1560 




K 


= 0.1554 




Qmi 


ti/a 


tf/a 


xVd.o.f. 


Gmi 


ti/a 


tf/a 


xVd.o.f. 


Gmi 


ti/a 


tf/a 


xVd.o.f. 


2.25(6) 


8 


12 


0.1 


2.35(5) 


8 


12 


0.2 


2.98(6) 


10 


16 


1.9 


2.25(7)* 


8 


16 


0.2 


2.36(5)* 


8 


16 


0.5 


3.00(7) 


10 


20 


1.6 


2.28(7) 


10 


12 


0.1 


2.37(6) 


10 


12 


0.002 


3.07(8) 


12 


16 


0.4 


2.26(9) 


10 


16 


0.2 


2.39(7) 


10 


16 


0.9 


3.09(9)* 


12 


20 


0.5 




ti/a 


tf/a 


xVd.o.f. 


Gb2 


ti/a 


tf/a 


xVd.o.f. 


Ge2 


ti/a 


tf/a 


xVd.o.f. 


0.036(10) 


8 


12 


0.3 


0.004(7) 


8 


12 


2 


0.049(15)* 


10 


16 


0.5 


0.038(11)* 


8 


16 


0.4 


0.003(7) 


8 


16 


1.4 


0.046(16) 


10 


20 


1.1 


0.037(15) 


10 


12 


0.6 


-0.010(11) 


10 


12 


1.0 


0.047(23) 


12 


16 


0.8 


0.042(17) 


10 


16 


0.5 


-0.012(13)* 


10 


16 


0.5 


0.038(26) 


12 


20 


1.2 


Rem 


ti/a 


tf/a 


xVd.o.f. 


Rem 


ti/a 


tf/a 


xVd.o.f. 


Rem 


ti/a 


tf/a 


xVd.o.f. 


-1.6(5) 


8 


12 


0.3 


-0.12(29) 


8 


12 


2.0 


-1.6(5)* 


10 


16 


0.6 


-1.7(5)* 


8 


16 


0.4 


-0.13(32) 


8 


16 


1.4 


-1.6(5) 


10 


18 


1.0 


-1.8(6) 


9 


16 


0.4 


0.01(41) 


9 


16 


1.0 


-1.5(5) 


10 


20 


1.2 


-1.6(6) 


10 


12 


0.6 


0.44(48) 


10 


12 


1.0 


-1.5(8) 


12 


16 


0.8 


-i.8(7) 


iO 


i() 


{).:> 


0.52(5(.))' 


iO 


i() 


().:> 


-i.2(8) 


12 


21) 






ti/a 


tf/a 


xVd.o.f. 


GC2 


ti/a 


tf/a 


xVd.o.f. 


Gc2 


ti/a 


tf 


xVd.o.f. 


















0.18(7) 


10 


16 


0.3 


















0.12(13)* 


12 


16 


0.1 


K = 0.155 


K 


= 0.1570 




K = 0.1562 


Gmi 


ti/a 


tf/a 


xVd.o.f. 


Gmi 


ti/a 


tf/a 


xVd.o.f. 


Gmi 


ti/a 


tf/a 


xVd.o.f. 


1.94(8) 


8 


12 


0.4 


1.92(5) 


8 


12 


0.2 


2.71(9) 


10 


16 


1.1 


1.92(9) 


8 


16 


0.9 


1.91(7)* 


8 


16 


0.1 


2.72(9) 


10 


20 


1.2 


1.95(10)* 


9 


15 


0.8 


1.93(8) 


9 


16 


0.3 


2.82(13) 


12 


16 


0.2 


2.01(12) 


10 


12 


0.02 


1.95(9) 


10 


12 


0.2 


2.83(13)* 


12 


18 


0.3 


1.94(13) 


10 


16 


1.1 


1.94(10) 


10 


16 


0.1 


2.84(14) 


12 


20 


0.6 


Ge2 


ti/a 


tf/a 


xVd.o.f. 


Ge2 


ti/a 


tf/a 


xVd.o.f. 


Ge2 


ti/a 


tf/a 


xVd.o.f. 


0.074(24) 


8 


12 


0.3 


0.022(9 ) 


8 


12 


0.2 


0.077(24)* 


10 


16 


0.5 


0.072(25) 


8 


16 


0.3 


0.023(10) 


8 


16 


0.4 


0.075(25) 


10 


20 


0.8 


0.078(33)* 


9 


15 


0.3 


0.023(13)* 


9 


15 


0.2 


0.105(47) 


12 


16 


0.2 


0.077(37) 


10 


12 


0.2 


0.027(14) 


10 


12 


0.1 


0.097(50) 


12 


18 


0.6 


0.071(38) 


10 


16 


0.3 


0.029(16) 


10 


16 


0.4 


0.093(50) 


12 


20 


0.8 


Rem 


ti/a 


tf/a 


xVd.o.f. 


Rem 


ti/a 


tf/a 


xVd.o.f. 


Rem 


ti/a 


tf/a 


xVd.o.f. 


-3.9(1.3) 


8 


12 


0.2 


-1.1(5) 


8 


12 


0.2 


-2.9(9)* 


10 


16 


0.4 


-3.8(1.3)* 


8 


16 


0.2 


-1.2(5) 


8 


16 


0.4 


-2.8(10) 


10 


20 


0.8 


-4.0(1.7) 


9 


16 


0.2 


-1.2(7)* 


9 


16 


0.4 


-3.7(1.7) 


12 


16 


0.3 


-3.8(1.9) 


10 


12 


0.2 


-1.4(7) 


10 


12 


0.2 


-3.4(1.7) 


12 


18 


0.6 


-3.6(1.9) 


10 


16 


0.2 


-1.5(9) 


10 


16 


0.4 


-3.2(1.8) 


12 


20 


0.9 


GC2 


ti/a 


tf/a 


xVd.o.f. 


GC2 


ti/a 


tf/a 


xVd.o.f. 


GC2 


ti/a 


tf/a 


xVd.o.f. 










0.11(3) 


8 


12 


0.2 


0.25(9) 


10 


15 


0.2 










0.11(3)* 


8 


15 


0.3 


















0.12(5) 


10 


12 


0.4 


0.17(21)* 


12 


15 


0.1 










0.13(5) 


10 


15 


0.3 
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TABLE III: Results for Qmi, Qe2 in units of natural magnetons (e/2mjv) and for the ratio —Qe2/Qmi- 









Quenched: l3 = 


= 6.0, 16^ X 32, q = (27r/16a, 0, 0) 






K 


p at rest 


A at rest 


p at rest 


A at rest 


p at rest 


A at rest 


A at rest 




Ma) 


nW 


Ma) 




Ma) 


Ma) 
f B2 


yE2 


Rem (%) 


yc2 


0.153 


2.37(7) 


2.52(7) 


2.25(7) 


2.26(7) 


0.046(13) 


0.038(11) 


0.035(12) 


-1.8(8) 


-1.7(5) 




0.154 


2.24(8) 


2.41(9) 


2.12(8) 


2.14(8) 


0.054(19) 


0.054(19) 


0.045(18) 


-2.4(9) 


-2.3(7) 




0.155 


2.15(14) 


2.27(11) 


1.95(10) 


1.99(10) 


0.068(28) 


0.078(33) 


0.071(24) 


-3.3(1.4) 


-3.8(1.3) 










Quenched: /3 = 


-- 6.0, 32^ X 64, q = (27r/16a, 0, 0) 






K 


p at rest 


A at rest 


p at rest 


A at rest 


p at rest 


A at rest 


A at rest 




Ma) 




Ma) 




Ma) 


Ma) 


yE2 


Rem (%) 


yc2 


0.1550 


1.87(10) 


1.99(10) 


1.78(7) 


1.76(7) 


0.048(20) 


0.049(40) 


0.056(41) 


-2.6(7) 


-0.95(48) 










Quenched: l3 = 


= 6.0, 32^ X 64, q = (27r/32a, 0, 0) 






K 


p at rest 


A at rest 


p at rest 


A at rest 


p at rest 


A at rest 


A at rest 




Ma) 




Ma) 




Ma) 
^£2 


Ma) 
f B2 


G(b) 

yE2 


Rem (%) 


Ma) 

yc2 


0.1554 


3.24(10) 


3.44(10) 


3.09(9) 


3.12(9) 


0.073(30) 


0.049(15) 


0.049(16) 


-2.2(9) 


-1.6(5) 


0.12(13) 


0.1558 


3.11(11) 


3.35(10) 


2.97(11) 


3.00(11) 


0.079(23) 


0.059(18) 


0.057(18) 


-3.0(1.2) 


-2.1(6) 


0.12(15) 


0.1562 


2.96(13) 


3.22(12) 


2.83(13) 


2.84(14) 


0.131(53) 


0.077(24) 


0.075(24) 


-4.1(2.0) 


-2.9(9) 


0.17(21) 








Nf 


= 2: p = 


5.6, 16^ X 32 q= (27r/16a, 0, 0) 








K 


p at rest 


A at rest 


p at rest 


A at rest 


p at rest 


A at rest 


A at rest 




Ma) 


M") 


Ma) 


Mb) 


Ma) 

yE2 


Ma) 

yE2 


p(b) 

yE2 


Rem (%) 


Ma) 

yc2 


0.1560 


2.51(6) 


2.65(6) 


2.36(5) 


2.36(5) 


0.031(8) 


-0.012(13) 


-0.011(12) 


-1.1(4) 


0.52(56) 




0.1565 


2.34(8) 


2.52(7) 


2.13(6) 


2.15(6) 


0.035(10) 


0.020(11) 


0.025(9) 


-1.3(5) 


-0.7(5) 


0.11(3) 


0.1570 


2.03(8) 


2.18(9) 


1.91(7) 


1.90(7) 


0.033(12) 


0.023(13) 


0.027(10) 


-1.3(7) 


-1.2(7) 


0.11(3) 


0.1575 


2.01(7) 


2.29(8) 


1.73(6) 


1.72(7) 


0.070(17) 


0.021(16) 


0.025(16) 


-3.5(9) 


-1.2(9) 


0.08(3) 



To extract a reliable value for the form factors we iden- 
tify the best plateau region having the largest possible 
range and ensure that, changing the fit range within this 
plateau region, produces results that remain within sta- 
tistical errors of each other. In TablelTllwe give the varia- 
tion of the mean values from changing the fit range of the 
plateaus for two representative k values for each lattice 
in the rest frame of the A. The quality of the plateaus 
and best fit ranges for the case when the nucleon is at 
rest is displayed in Figs.llUlandllllfor Qmi and in Fig. 1121 
for Ge2- In Fig. E| we display the plateaus for the large 
quenched lattice for Ge2 in the rest frame of the A in 
order to have one case of direct comparison to the num- 
bers chosen for the best plateau in Table ^ In Fig. H"!] 
we show the plateaus and best fit ranges for Rem in the 
rest frame of the nucleon. From the plateaus values given 
in Table ITU it can be seen that the dependence of the re- 
sults on the fit ranges are, in all cases, well within the 
statistical errors obtained from jackknife analysis. 

In order to access finite volume effects we calculate the 
transition form factors in the quenched theory on two lat- 



tices of size 16'^ x32 and 32'^ x 64 at the same momentum 
transfer. We show the results for Gmi, Qe2 and the ratio 
Ge2/Gmi in Fig. 121 The data for Gmi obtained on the 
large volume lie systematically below the ones obtained 
on the small lattice. The results for Ge2 as well as for 
the ratio Ge2/Gmi, have overlapping errors. On the two 
volumes, the difference in the plateau values is smallest 
in the ratio since volume effects tend to cancel in numer- 
ator and denominator. The plateau values for Gmi given 
in Table UTTl are systematically lower by 10-15% on the 
larger lattice. Therefore there is a volume correction on 
the results obtained on the small lattices of the order of 
10% which is comparable to the statistical error. The re- 
sults obtained on the larger lattice have negligible volume 
dependence if the volume correction scales inversely pro- 
portional to the volume. For Ge2 and the ratio Ge/Gmi 
which carry larger statistical errors the volume depen- 
dence is harder to access since the difference in the values 
for the two volumes is well within the statistical errors. 



The values for Gmi extracted from the fits to the 
plateaus are collected in Table Hill The overall quality of 



K value can be seen in Fig. ^1 for q = (27r/16a, 0, 0) and 
in Fig. ^2 for q — (27r/32a, 0, 0) for both the quenched 
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FIG. 10; Gj^ii at momentum transfer q = (27r/16a, 0, 0) in 
the rest frame of the nucleon. The three upper graphs on 
the left show the quenched resuhs for the lattice of size 16'^ 
at K = 0.153, 0.154 and 0.155. The lowest graph on the left 
shows quenched results for the lattice of size 32^. On the right 
we show the unquenched results at k = 0.1560, 0.1565, 0.1570 
and 0.1575. The photon couples to the quark at ti/a — 6 for 
the small lattices and at ti/a — 8 for the large lattice as shown 
by the arrow. The dashed lines show the fit range and bounds 
on the plateau value obtained by jackknife analysis. 



the lattice improves the identification of the plateau re- 
gion yielding reliable results for Gmi at the lighter quark 
masses. From the values given in Table Hill we conclude 



that, when the A is produced at rest, definitions Gm\ or 
^Mi yield the same values. For the case where the A 
carries momentum q along the x-direction which 
is extracted from matrix elements for which the A vec- 
tor index is in the x-direction, yields consistently larger 
values than those obtained from Q'}^\ ■ This difference in 

the values of Gm\ and Gmi volume independent and 
increases from about 6% for the heavier quarks to about 
15% for the lighter ones. It is also independent of the 
value of the momentum carried by the A. It would be 
interesting to allow momentum transfers in all directions 
to check if this difference will be reduced and also use an 
0(a) improved Dirac operator to check whether it is due 
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FIG. 11: Gmi ^^o'' the quenched theory on a lattice of size 
32=* X 64 at K = 0.1554, 0.1558 and 0.1562 with the nucleon 
at rest and for momentum transfer q — 27r/32a. The photon 
couples to the quark at fi/a = 8 as indicated by the arrow. 
The dashed lines show the fit range and bounds on the plateau 
value obtained by jackknife analysis. 



To look for sea quark effects on the value of Gmi we 
compare the quenched results on the small lattice to the 
unquenched results using the SESAM configurations at 
similar ratio of the pion mass to the rho mass. From the 
values of pion to rho mass ratios given in TablelHone thus 
compares quenched results at k = 0.153, 0.154 and 0.155 
to the unquenched results at k = 0.1560, 0.1570 and 
0.1575, respectively. Unquenching leads to a stronger 
quark mass dependence increasing the value of Gmi at 
the heaviest quark mass and reducing it at the two lighter 
quark masses. 

In Fig. 1121 we show the quenched and unquenched re- 
sults for g'e2 at q = (27r/16a, 0, 0) in the lab frame of the 
nucleon. The plateau region is limited to a a few time 
slices for the quenched lattice whereas for the SESAM 
lattice for twice the statistics the fits can be extended 
over a larger time interval. For a static A both G^e2 ^^'^ 
^7^2 show a similar behaviour and the best plateau fit 

ranges are given, for 5^2, in Table HTl for representative 
K values. As we already mentioned, although the central 
value for Ge2 decreases by going to the larger lattice, 
the observed decrease is well within our statistics. Sim- 
ilarly unquenching systematically reduces the value of 
Ge2- However with our statistics this reduction remains 
within errors, for the three k values that correspond to 
similar ratios of pion to rho mass. Plateau identification 
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display results a,t k — 0.155 for both the large and small 
quenched lattice at equal momentum transfer. This im- 
provement is also seen in Fig. 1131 where we show results 
for the large quenched lattice at q = (27r/32a, 0, 0), even 
though lighter quarks are used. Here part of the noise 
reduction is due to having a smaller value of the momen- 
tum transfer. This enables to check the stability of our 
fits by changing the fit range as given in Table Hll 
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FIG. 12: g'j^^. The notation is as in Fig. ITHl 



Our lattice results for the ratio —G'e'z/S'mi ^^i^" 
played for the quenched and the unquenched theory in 
Fig. El in the rest frame of the nucleon. The general 
trend is that this ratio becomes more negative as we ap- 
proach the chiral limit. The plateau values seen in the 
figure and given in Table IIIII give no indication of an 
increase in this ratio as we unquench. This may mean 
that pion cloud contributions, expected to drive this ra- 
tio more negative, are suppressed at these large quark 
masses. However one must keep in mind that part of the 
pion cloud is taken into account in the quenched theory 
since, using relativistic quarks includes backward propa- 
gation in time and thus effectively pionic contributions. 
A study with lighter pions is required in order to assess 
the importance of pion contributions due to sea quarks on 
the value of this ratio and thus on nucleon deformation. 

Finally in Fig. we show the results for CJ^2 with A 
static for the large quenched lattice for which the ratio i?^ 
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to allow plateau identification. Although an evaluation 
of 

Gc2 is ^iso carried out for the SESAM configurations 
a large enough plateau range could not be identified for 
all K values. Poor plateau identification leading to a lim- 
ited fit range may introduce systematic errors that are 
not included in the errors of the extracted values given 
in Table IhH However for the cases where we can fit over 
4 time slices far enough from the time of the current in- 
sertion the mean value of positive at all k values. 
Therefore even though the systematic errors are largest 
for this form factor, which means that the jackknife er- 
rors given in Table Hill provide an underestimation of the 
actual error, our lattice results favour a negative value 
for the ratio —Gc2/Gmi in agreement with experiment. 



IV. CHIRAL EXTRAPOLATIONS 

In order to obtain physical results we need to extrapo- 
late the lattice data to the chiral limit. Chiral perturba- 
tion theory has been applied to calculate the transition 
form factors |l3j but the range of validity is limited to 
very small quark masses and very low momentum trans- 
fers and therefore it cannot be used in the current analy- 
sis 25|. Since the nucleon or the A carry a finite momen- 
tum we expect chiral logs that appear at next-to-leading 
order in chiral perturbation theory to be suppressed for 
the momentum transfers studied in this work. There- 
fore a behaviour proportional to the pion mass squared 
should be appropriate all the way to very near the chi- 
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FIG. 14: The ratio 



(a) 



in the lab frame of the nu- 



cleoli. The three upper graphs on the left show results 
for the quenched theory on a lattice of size 16^ X 32 for 
q = (27r/16a,0,0) at k ^ 0.153, 0.154 0.155 and the lowest 
shows results for a lattice of size 32^ X 16 for the same value of 
q at K = 0.155. On the right we show results for the SESAM 
configurations at k = 0.1560, 0.1565, 0.1570 and 0.1575 for 
q = (27r/16a, 0, 0). The current couples to the quark at time 
ti/a = 6 for the lattices of temporal size 32a and at ti/o = 8 
for the lattice of temporal size 64a as shown by the arrow. 



pion mass, ttIt^, squared is seen from Figs. 1161 1171 and 
IT51 The results shown in these figures are for the kine- 
matical case where the A is stationary and remains 
almost unchanged as the pion mass approaches the chi- 
ral hmit. The errors on the lattice data are obtained by 
jackknife analysis and are purely statistical. Given the 
fact that we have results only at three values of k in the 
quenched case and four in the unquenched case a linear 
fit is the best option we have. For the quenched calcu- 
lation the error on the extrapolated value at the chiral 
limit that we quote in Table IIVI is obtained by doing a 
complete jackknife analysis for the three n values. In the 
unquenched case, configurations at each quark mass are 
obtained from different Monte Carlo simulations, and a 
standard ^^-minimization procedure, which assumes un- 
correlated data points at each quark mass, can be ap- 
plied. 

Extrapolating linearly in mi we obtain for both 
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FIG. 16: Top: G'^^'jI, middle: S^'t] and bottom:EMR in % 
versus in lattice units in the rest frame of the A for the 
quenched 16'^ x 32 lattice. The dashed line is fit to a + 6m^ 
and the dashed-dotted to a constant with the chiral value ob- 
tained shown with a filled square and star for the two Ansatze 
respectively. 



15 




0.00 0.05 



0.10 0.15 

(m a)2 



0.20 0.25 



^ 3.2 
^ 3.0 
\ 2.8 
^ 2.6 
2.4 



^ -12 
.10 

CN 

\ .08 

(D 

.06 

(N 
UJ 

U .04 

-1 

-2 
-3 
-4 
-5 



32^x64 N =0 



0.00 



0.02 0.04 0.06 
\2 



0.08 



(m a)^ 



FIG. 17: The same as in Fig. 1161 but for the unquenched 
theory. 



FIG. 18: The same as in Fig.[THlbut for the quenched 32^ x 64 
lattice. 



which is the most accurately determined quantity. All the 
data points fall nicely on a straight line even in the un- 
quenched case where the data at different quark masses 
are uncorrelated since they are obtained using a differ- 
ent set of configurations. A linear Ansatz also fits well 
the quenched results for Qe2 and Rem as can be seen in 
Figs. El El and El For the corresponding unquenched 
results, although a linear Ansatz still provides a good fit 
giving a x^/d.o.f. = 1.1 for Qei and x^/d.o.f. = 0.6 for 
Rem, the unquenched results for Qe2 and Rem show 
very weak mass dependence at the three lightest quark 
masses. Chiral perturbation theory suggests a similar 
mass dependence for Qe2 as that for Qmi- To check this 
mass dependence an unquenched calculation with higher 
statistics on a larger lattice to avoid finite size effects 
that, especially at the lightest quark mass (k — 0.1575), 
can be significant is called for. Given this weak mass de- 
pendence we also fit Qe2 and Rem to a constant which 
gives a lower value for Qe2 and the absolute magnitude 
of Rem at the chiral limit. For comparison we include 
also for the quenched results the value obtained using a 
constant fit even though a linear fit in is favored by 
the data. Although the chiral absolute value of Rem is 
reduced by modifying the fitting Ansatz to a constant. 
Rem remains negative in all cases. 

For kinematics where the nucleon is at rest changes 
as we approach the chiral limit. This change is partic- 
ularly severe for the large quenched lattice where one 
would need the extrapolation of form factors computed 
at that decrease by 50%. In phenomenological studies 
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FIG. 19: The EMR ratio in % in the rest frame of the nu- 
cleon. Top for the quenched 16'^ x 32 lattice, middle for the 
unquenched theory and bottom for the quenched 32'^ x 64 lat- 
tice . The dashed line is fit to a -I- 6m^ and the dashed-dotted 
to a constant with the chiral value obtained shown with a 
filled square and star for the two Ansatze respectively. 
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linear fit yielding results that are consistent with each 
other. We show only results obtained in the rest frame 
of the A for the large quenched lattice where we have 
the best signal. For this lattice and kinematics at = 
0.13 GeV^ the values we extract from the linear Ansatz 
are q'^^I = 0.20(32) and Q%\ = 0.48(34) in units of natu- 
ral magnetons (e/2mAr). Changing the fitting range can 
lead to changes of the mean value of about 30% at each 
K value. Although the mean values of ^"^^ system- 
atically larger than those of ^ values they are 
always within statistical errors. Despite the large statis- 
tical and systematic uncertainties on this quantity the 
data favour more a negative than a positive CMR ratio. 
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a + brn^ and the dashed-dotted to a constant with the chiral 
value obtained shown with a filled square and star for the two 
Ansatze respectively. 



by 



(33) 



where Ra{Q^) for a = Ml, E2 and C2 measures the de- 
viations from the proton electric form factor G^{Q^) — 
1/(1 -I- (3^/0.71)^. Usually experimental data are fit- 
ted by taking RmAQ^) = RE2iQ^) = Rc2{Q^) = 
1 -I- aexp(— (3^7) |23|. There are not enough data to si- 
multaneously fit the and the quark mass dependence. 
We are in the process of studying the dependence of 
these form factors using the fixed source sequential tech- 
nique [2^ |2^. However, assuming that the above phe- 
nomenological Ansatze provide a good description for the 
dependence of Gmi and Ge2, we expect the de- 
pendence to cancel in the ratio Rem- Therefore a chiral 
extrapolation can be performed for Rem also in the case 
of the nucleon being at rest. A constant and a linear fit 
in toJ are shown for the three lattices in Fig. ^| Once 
more the constant provides a lower limit for the absolute 
magnitude of Rem and still leads to a negative value in 
the chiral limit in all cases. 

Finally we can not observe any quark mass depen- 

»C2 / ^Ml 



dence for Q'f^l and QY^L as well as for the ratio 



,(6) 



^C2 



^C2 



given the present statistical uncertainties as can be seen 



V. DISCUSSION 

In Table IIVI we list the chiral extrapolated values of 
the Ml and E2 form factors obtained in the rest frame 
of the A as well as of the ratio Rem obtained for both 
kinematics. The values given are the result of using a lin- 
ear Ansatz in rn^ to extrapolate to the chiral limit and 
the quoted errors are only statistical. As discussed in 
section III and IV there is a number of systematic errors 
that we must consider when comparing the results of Ta- 
ble llVl to the experimental ones. We summarize here the 
source of systematic errors on the form factors and their 
ratios. Finite volume effects are estimated by performing 
a quenched calculation at the same value of momentum 
transfer on a lattice of size 16'^ x 32 and 32^ x 64 at 
(3 — 6.0. Assuming a 1/volume dependence for Qmi the 
results show a 10-15% correction for the data obtained 
on the small lattice and negligible for the larger lattice. 
The statistical errors on Qe2 are too large to enable any 
volume correction to be extracted. However one would 
expect a similar volume dependence as for Qmi- On the 
other hand the volume dependence largely cancels in the 
ratio and so we expect both EMR and CMR ratios to 
have negligible volume dependence. The variation in the 
mean values obtained by varying the plateau ranges is 
smaller as compared to the statistical errors as can be 
seen from Table ^1 where we have given the values ob- 
tained using different fit ranges for representative k val- 
ues for all the lattices studied in this work. We look 
for unqucnching effects by comparing quenched and un- 
qucnchcd results on lattices of similar physical volume 
and lattice spacing as well as pion to rho mass ratio. For 
pion masses in the range of 800-500 MeV unquenching 
effects are within statistical errors. Extrapolation of the 
lattice results to the chiral limit represents the biggest 
uncertainty that can only be eliminated by evaluating 
the form factors closer to the chiral limit. Avoiding any 
chiral extrapolations the results on the two quenched and 
on the unqucnchcd lattices at the lightest quark mass give 
a ratio Rem in the range of about (-1 to -5)%, which is 
in accord with the range obtained in experimental mea- 
surements. Both quenched and unquenched calculations 
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of about 0.1 fm. Estimating finite lattice spacing effects 
would require repeating the calculation at two larger val- 
ues of (3 keeping the physical volume constant and at- 
tempting a continuum extrapolation. This is beyond our 
current computer resources and no systematic error due 
to the finite lattice spacing can be estimated. 

In the last column of Table IIVI we give the results 
for the ratio Rem obtained in the SL model using an 
effective Hamiltonian defined in the subspace of ttN 
7A^ © A |2]|. Within the SL model Rem is calculated 
in the rest frame of the A as a function of for bare 
and dressed vertices. It is interesting that for bare ver- 
tices they obtain non-zero values for Rem pointing to 
the fact that a non-zero electric quadrupole amplitude 
implying nucleon deformation is obtained even without 
pions. This is in agreement with lattice results on the 
rho deformation which is observed in the quenched ap- 
proximation without pion contributions 5]. Comparing 
the value obtained in the SL model at = 0.64 GeV^ 
and 0.13 GeV^ using bare vertices to our quenched lattice 
result via Eq.El 2i| we see that quenched QCD produces 
a more negative value for Rem than that obtained in 
the SL model. One must however take into account that, 
in the baryonic sector, quenching still includes part of 
the pion cloud due to quarks propagating backwards in 
time. These are absent in the SL model with bare ver- 
tices. On the other hand unqucnched lattice results at 
quark masses that correspond to pion masses in the range 
500-800 MeV do not influence the value of Rem- This 
implies that pion cloud contributions from sea quarks at 
these pion masses are not large. It remains an open issue 
whether this means that pionic contributions from back- 
going quarks are also small, in which case the difference 
between quenched lattice and SL model results cannot be 
explained solely by pion contributions from back-going 
quarks. Results for Rem obtained in the SL model with 
fully dressed vertices are twice as negative as those ob- 
tained with bare vertices showing, within this model, the 
importance of pionic contributions to deformation. Our 
unqucnched results, obtained here for rather heavy pions, 
show no significant enhancement which may reflect that 
qq creation is suppressed. Since it is for light pions that 
large pionic contributions are expected it is imperative 
to repeat the calculation with dynamical quarks closer 
to the chiral limit to study sea quark contributions. 

The values of the ratio Rem given in Table IIVI can 
be compared to those measured at various values of 
in recent experimental searches for nucleon deformation. 
The lattice results for Rem are plotted in Fig. [Altogether 
with the recently measured experimental data. In the 
figure we also included the Particle Data group value 
for Rem at = Once more what can be seen 
from this figure is that quenched and unqucnched results 
are within errors. The unqucnched results at our two 
available momentum transfers are in agreement with the 
experimental results. We stress that all the errors shown 
on this figure are statistical. For the experimental results 



TABLE IV: Results for Qmi, Ge2 in units of e/2mjv and for 
the ratio Rem = —Qe2/Qmi in % extrapolated to the chi- 
ral limit. All errors given are statistical. The last column 
gives t he p rediction for Rem in % obtained within the SL 
model |21ll where for the quenched case we quote their val- 
ues without pion cloud contributions and for the unqucnched 
theory we give their fully dressed results for the case where 
the A is produced at rest. In the last two columns we give 
the experimental results at similar values of . 
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FIG. 21: The ratio Rem in % versus at the chiral 
limit. Filled triangles are quenched and open triangles are 
unquenched results. The filled square is the Particle Data 
Group result at Q'^ = [2^. the crosses are from ref. 2] and 
the star from ref. The errors shown on the experimen- 
tal results are statistical. The lattice values shown are those 
obtained using a linear fit in to extrapolate to the chiral 
limit. In the unquenched case the error bars are increased so 
that an extrapolation to the chiral limit using a constant lies 
within the shown error bands. 



error. For the lattice results the most severe systematic 
error comes from assuming a linear Ansatz in for the 
extrapolation to the chiral limit. 

From the lattice data obtained in this work Rsm 
can only be estimated in the chiral limit for the large 
quenched lattice. Extrapolating the results obtained 
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TABLE V: /mi, /b2, A1/2 and A3/2 in units of GeV"^''^ 



the values at = 0, 



-0.135 ± 0.006 GeV"i/2 



lattice 


GeV^ 


fMl 


/b2 


Al/2 


A3/2 


16^ X 32 


0.64 


0.266(22) 


-0.016(7) 


-0.109(13) 


-0.244(20) 


32^ X 64 


0.13 


0.296(23) 


-0.012(4) 


-0.129(12) 


-0.267(20) 


SESAM 


0.53 


0.194(15) 


-0.0078(36) 


-0.085(8) 


-0.175(13) 



and A 



3/2 



-0.255 ± 0.008 GeV-1/2, as given by the 



Particle Data Group. 



VI. CONCLUSIONS 



Q2 = 0.13 GeV2 we find -g^c2/SMi = (-7.7 ± 11)% 
which leads, via Eq. [HI to Rsm — (—1-2 ± 1.7) % where 
only the statistical error is given. At = 0.126 GeV^ 
the experimental value of Rsm is (—6.5 ±0.2 ±2.5) % 0. 
As already mentioned, the lattice evaluation of this quan- 
tity is affected by large systematic and statistical errors, 
which must be studied before a more accurate determi- 
nation can be obtained. 

Finally, in order to facilitate comparison with exper- 
iment and phenomenology, not only for the ratios but 
also for Gmi and Qe2 separately, we give the relationship 
between the Sachs form factors studied in this work and 
the electromagnetic transition amplitudes /mi and fE2 
in the rest frame of the A 1211: 



(Ml 



2m 



N 



1/2 



1 - q^/{mpf + rriAy 



1/2 



(34) 



/b2 — — 



q|nT-ANi/2 



2mN niN 



Ge2 

(35) 

with e = ■\/47r/137. Since these are continuum relation- 
ships to obtain results for /mi and fE2 from the values 
given in Table IIVI one uses the physical nucleon and A 
mass. We give the values for Jmi and fE2 in Table IVl 
They are are related to the helicity amplitudes A1/2 and 
A3/2 by 



/a/1 = -2(^1/2 + v^Ag/a) 
/s2 = -2(^1/2 - —^^3/2) 



(36) 



As can be seen from Table HVl the values of the helicity 
amplitudes for closest to zero are in agreement with 



The matrix element for the transition jN — s- A is 
computed in lattice QCD both for the quenched the- 
ory and for two dynamical Wilson fermions in a first 
attempt to study unquenching effects on these form fac- 
tors. The dominant magnetic dipole form factor is calcu- 
lated with statistical accuracy of about 10%. The electric 
quadrupole is suppressed by an order of magnitude and 
it is calculated to a statistical accuracy of about 50%. 
Unquenching tends to decrease the values of Gmi and 
gE2- The ratio of these form factors provides a direct 
comparison to the experimentally measured ratio Rem- 
We find a negative value for Rem of the order of a few 
% in accord with experiment. For pions of mass in the 
range of about 800 — 500 MeV we obtain no evidence 
for an increase in the value of Rem as we unquench. It 
is expected that pion cloud contributions are suppressed 
for these heavy quarks and therefore it is important, in 
future studies, to use lighter dynamical quarks for the 
evaluation of these form factors. Large statistical and 
systematic errors prevent at this stage a determination 
of the Coulomb quadrupole form factor. A detailed study 
of lattice artifacts will be needed for better control of sys- 
tematic errors, before a more accurate extraction of the 
Coulomb form factor and of the ratio Rsm is possible. 
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